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ABSTRACT 



We investigate the radio afterglows of gamma-ray bursts (GRBs) and hyper- 



' novae (HNe) at high redshifts and quantify their detectabihty, as well as their 

^ ■ potential usefulness for 21 cm absorption line studies of the intergalactic medium 

^ ■ (IGM) and intervening structures. We examine several sets of source and envi- 

■ ronment model parameters that are physically plausible at high redshifts. The 

p i' radio afterglows of GRBs would be detectable out to 2; ~ 30, while the energetic 

Q ■ HNe could be detectable out to z ~ 20 even by the current Very Large Array 

^ • (VLA). We find that the 21 cm absorption line due to the diffuse neutral IGM is 

^ ■ difficult to detect even by the proposed Square Kilometer Array (SKA), except 
for highly energetic sources. We also find that the 21 cm line due to collapsed 

^ . gas clouds with high optical depth may be detected on rare occasions. 

Subject headings: gamma rays: bursts — intergalactic medium — radio lines: 
general 



1. Introduction 

Gamma- Ray Bursts (GRBs) are potentially powerful probes of the early universe (e.g., 
Miralda-Escude 1998; Meszaros& Rees 2003; Barkana & Loeb 2004). The substantial evi- 
dence that massive stellar collapses lead to long duration GRBs suggests that these outbursts 
probably occur up to the highest redshifts of the first generation of stars (c.f., Meszaros 2002; 
Zhang & Meszaros 2003). Their high luminosities make them detectable in principle out to 
redshifts z ~ 100 (Lamb & Reichart 2000), while their afterglows may be observable up 
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to z ~ 30 in the infrared, X-ray (Ciardi & Loeb 2000; Gou, Meszaros, Abel, & Zhang 
2004) and radio bands (loka 2003; Inoue 2004). About 25% of all GRBs detected by the 
upcoming Swift satellite are expected to be at 2; > 5 (Bromm & Loeb 2002). Such high 
redshift GRBs may have already been detected by the BATSE satellite, as suggested by the 
empirical distance indicators (Fenimore & Ramirez-Ruiz 2000; Norris, Marani, & Bonnell 
2000; Murakami et al. 2003; Yonetoku et al. 2004) and their theoretical interpretation 
(e.g., loka & Nakamura 2001; Yamazaki, loka, & Nakamura 2004). The first generation 
of stars are expected to be very massive (Abel, Bryan, & Norman 2002; Bromm, Coppi, & 
Larson 2002; Omukai & Palla 2003), so that their outcome might be GRBs which are an 
order of magnitude more luminous than their low redshift counterparts (e.g.. Fryer, Woosley, 
& Heger 2001; Heger et al. 2003). 

One of the major uncertainties in modern cosmology is the reionization history of the 
universe, which has significant implications for the formation of the first stars and their 
feedback on the subsequent evolution (e.g., Barkana & Loeb 2001; Miralda-Escudc 2003). 
The analysis of Lya spectra in the highest redshift quasars indicate that the reionization 
has essentially reached completion at around z ~ 6 (e.g.. Fan et al. 2002a), while the 
WMAP polarization data imply an onset of the reionization process at much higher redshifts 
2; ~ 17 ± 5 (Kogut et al. 2003; Spergel et al. 2003). Although other methods, such as the 
temperature of the Lya forest (Theuns et al. 2002; Hui & Haiman 2003), the Hn region 
size around quasars (Wyithe & Loeb 2004; Mesinger & Haiman 2004), and the luminosity 
function of the Lya emitters (Stern et al. 2004; Malhotra & Rhoads 2004) also constrain the 
reionization history, it is still unclear even whether the reionization is a continuous process 
or whether it occurred in two episodes (Gen 2003; Wyithe & Loeb 2003). 

21 cm tomography would be one of the promising ways to reveal the reionization history 
of the universe^ (Tozzi et al. 2000; Shaver et al. 1999), by mapping the neutral hydrogen in 
the spatial and redshift space via the redshifted 21 cm hue (Madau, Meiksin, & Rees 1997; 
Scott & Rees 1990; Hogan & Rees 1979). The 21 cm line is associated with the transition 
between the singlet and triplet hyperfine levels of the hydrogen atom at the ground state 
(e.g.. Field 1959a). The 21 cm line may be detectable in emission or absorption against the 
cosmic microwave background (GMB) (e.g., Giardi & Madau 2003; Furlanetto, Sokasian, & 
Hernquist 2004; Zaldarriaga, Furlanetto, & Hernquist 2004; Ihev et al. 2002). Alternatively, 
the 21 cm line may be also observed in absorption against luminous background sources 
(Garilli, Gnedin, & Owen 2002; Furlanetto & Loeb 2002). But this method requires 



^Other methods using GRBs may include the infrared spectroscopy (Miralda-Escude 1998; Barkana & 
Loeb 2004) or photometry (Inoue, Yamazaki, & Nakamura 2004) and the dispersion measure (loka 2003; 
Inoue 2004). 
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luminous sources of radio continuum at high redshifts, which remain largely speculative. 
The GRB radio afterglows would be a natural candidate. However their radio fluxes at high 
redshift have not been estimated in detail. Furlanetto & Loeb (2002) suggested that GRB 
afterglows are too dim, but they did not explicitly calculate the self-absorption frequencies, 
nor did they include the reverse shock emission. Another type of related energetic radio 
source which has not been investigated in this respect are hypernovae (HNe), which may be 
as frequent as GRBs. 

In this paper, we investigate the GRB and HN radio afterglows at high redshifts, in- 
cluding also the possibly of more energetic versions of these objects which may be associated 
with a population III of stars. We first estimate the maximum redshifts out to which these 
afterglows could be detected with the Very Large Array (VLA),^ Low Frequency Array (LO- 
FAR),^ and Square Kilometer Array (SKA).^ This can be important for broadband afterglow 
fits aimed at determining total energies and physical parameters (e.g., Panaitescu & Kumar 
2002) of bursts, out to very high redshifts. We also quantify the detectability of the 21 
cm absorption line in the afterglows. Since the afterglow model parameters, such as the 
energy and duration of the GRBs as well as the ambient density may be different at high 
redshifts, we examine several possibilities of the model parameters. However, even with the 
SKA, which has the best sensitivity among the proposed telescopes, and with our improved 
spectral-temporal evolution model, we find that the chances are low for the detection of 
21 cm absorption fines from GRBs and HNe, except in the case of exceptionally energetic 
sources. This quantification exercise should be useful for the assessment and planning of 
future observing projects with advanced facilities and detectors. 

The paper is organized as follows. In § 2, we calculate the GRB and HN radio afterglows 
at high redshift, based on the afterglow model summarized in § A. By comparing the 
afterglow fiuxes with the VLA, LOFAR and SKA sensitivities, we discuss the maximum 
redshift out to which the radio afterglow of these sources would be detectable. We calculate 
the 21 cm line absorption in the IGM and intervening structures at z > 6 in § 3 and at 
^ < 6 in § 4. Our conclusions are summarized in § 5. Throughout the paper we use a ACDM 
cosmology with {Q^, flf,, fib, h) = (0.27, 0.73, 0.044, 0.71) (Spergel et al. 2003). 



^http:/ /www. nrao.edu/ 
■^http:/ /www. lofar.org/ 
http : / / www . skatelescope . org/ 
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2. Radio fluxes of high redshift afterglows 

We model the GRB afterglow source as a relativistic shell expanding into a homogeneous 
interstellar medium (ISM) with particle number density n at redshift z. The shell initially 
has an isotropic equivalent energy E, a Lorentz factor 70, an opening half- angle 6 and a width 
in the source frame cT{l + z)~^, where we assume 7^^^ < 9 and the shell width being related 
to the observed GRB duration T (Kobayashi, Piran, & Sari 1997). The true energy is given 
by Ej — 9'^E/2. Two shocks are formed: a forward shock heating the ISM and a reverse 
shock decelerating the shell. At these shocks electrons arc accelerated and magnetic fields 
are amplified, leading to the synchrotron afterglow emission. We assume that accelerated 
electrons have a power-law distribution of the Lorentz factor 7e as iV(7e)(i7e oc 7~^(i7e for 
7e > 7m (Sari, Piran, & Narayan 1998). We also assume that fractions eg,/ {^e,r) and 
^Bj {^B,r) of the shock energy go into electrons and magnetic fields, respectively, at the 
forward (reverse) shock, where the subscripts / and r indicate the forward and reverse 
shock, respectively. We also assume adiabatic shocks. 

Under the above assumptions, we can calculate the spectra and light curves of the 
afterglows as summarized in § A. There are 10 model parameters: E, 6, n, p, Sgj, ^bj, 
^e,r, ^B,r, 7o, T / {1 + z) , but we assume Cej — ee,r for simplicity, so the actual number 
of parameters is 9. We consider the possible difference between the magnetic field in the 
forward and reverse shocks, by using TZb — {'^B,r/ ^bjY^"^ in equation (A24) instead of eB,r-i 
since the ejected shell may be endowed with magnetic fields from the central source (e.g., 
Zhang, Kobayashi, & Meszaros 2003). We also take the sideway expansion of the jet, the 
non- relativistic regime, and the reverse shock emission into account (see § A). 

In the following we examine several sets of model parameters that are physically moti- 
vated at high redshift, as summarized in Table 1. We estimate the maximum redshift out to 
which these afterglows could be detected with the VLA, LOFAR and SKA. We also evaluate 
the peak fluxes of the afterglows at the observed frequency around u ~ 100 MHz to discuss 
the detectability of the 21 cm absorption line later in § 3 and § 4. 

2.1. Standard GRB case 

First we adopt typical parameters inferred from the broadband afterglow fitting (e.g., 
Panaitescu & Kumar 2002; Zhang, Kobayashi, & Meszaros 2003): the isotropic equivalent 
energy E = 10^''' erg, the opening half-angle 6 = 0.1, the ISM density n = 0.1 cm^^, the 
spectral index p = 2.2, the initial Lorentz factor 70 = 200, the duration in the source frame 
T/(l + z) = 10 s, and the plasma parameters eg = 0.1, = 0.01, TZb = 1. 
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Figure 1 shows the total fluxes of the forward and reverse shock emission at the observed 

frequency z/ = 5 GHz and the observed times, t = 1 hr, 1 day, 10 days and 100 days, as a 
function of redshift z. The 5(t sensitivities of the VLA and SKA are also shown. Wc can see 
that the radio afterglows can be detected up to ^ ~ 30 even by the current VLA. Here we 
use the sensitivity. 
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assuming the signal-to-noise ratio (SNR) 5, the integration time tint = 1 day, the band width 
Aiy = 50 MHz, ^eff/T^ys ~ 2 x 10^ cm-^ K'^ for the VLA, and A^s/T^ys ~ 2 x 10^ cm-^ 
for the SKA. Note that LOFAR will has the frequency range from ~ 10 to ~ 250 MHz, 
and not at ~ 5 GHz. In Figure 1, the forward shock emission usually dominates the reverse 
shock one. We can also see that the redshift dependence is rather weak. 

Figure 2 shows the hght curves of the standard GRB afterglows at 2; = 6 and u ~ 200 

MHz and at z = 13 and u ~ 100 MHz, which correspond to the redshifted 21 cm band. The 
5a sensitivities of the VLA, LOFAR and SKA are also shown. We see that LOFAR and SKA 
can detect the standard afterglows at low frequencies ~ 100 MHz and high redshift z ~ 10, 
but it may be difficult to detect them with the VLA. Here we assumed an integration time 
of one-third of the observed time t, a band width Ai> = 50 MHz, Af^g/Tsys ~ 3 x 10^ cm~^ 
K-i for the VLA, ^e/Tsys ~ 4 x 10^ cm-^ K'^ for the LOFAR, and A^s/Tsy, ~ 5 x 10^ 
cm-2 K-^ for the SKA. Note that Aeg/Tsys at ~ 100 MHz is lower than that at ~ 5 GHz 
because of the galactic background. 

From Figure 2 one sees that the peak flux is about ~ 1-10 //Jy at the redshifted frequency 
of the 21 cm line, u — 1420/ (1 -|- z) MHz, for z >6. At these frequencies, the forward shock 
is brighter than the reverse shock, and the dependence on the redshift is rather weak. For 
most of the time the afterglow flux is below the self-absorption frequency so the flux 
is limited by the blackbody value in equation (A20). Since the flux is proportional to the 
electron temperature and the reverse shock has lower temperature than the forward shock, 
the reverse shock is dimmer than the forward shock for our parameters. The first break 
in the forward shock emission at t ~ 10~^ day corresponds to the shock crossing time tx 
in equation (A2), the second break at t ~ 10 day corresponds to the jet break time te 
in equation (A8), the third break at i ~ 10^ day is due to the characteristic synchrotron 
frequency i/^j crossing the observed frequency v, the final peak is due to the self-absorption 
frequency Uaj crossing the observed frequency 1/, and the non-relativistic phase begins soon 
after the final peak (see § A). The self-absorbed fiux depends on the density as F,^,bb oc 
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at tx < t < te ii v < v^j < Ucj from equation (A20). Thus the flux is higher for a lower 
density. However, if the density is too low, the synchrotron flux in equations (A14) and 
(A15) falls below the blackbody flux in equation (A20) and the observed flux decreases as 
the density decreases, since -Fj^^maxj oc n^/^ at tx < t < ii v < z/^j < Vcj- We find that 
the flux is maximized at n ~ 10~^ cm~^ for our parameters (see also § 2.3). 

2.2. Energetic GRB case 

The first generation of stars are expected to be metal-free and very massive (Abel, Bryan, 
& Norman 2002; Bromm, Coppi, & Larson 2002; Omukai & Palla 2003). If the collapse 
of some of these population III stars leads to the launch of relativistic jets, they could lead 
to an order of magnitude more energetic GRBs than their low redshift counterparts (e.g., 
Meszaros& Rees 2003; Fryer, Woosley, & Heger 2001; Heger et al. 2003). Here we study a 
hypothesized energetic GRB case with isotropic equivalent energy E — 10^^ erg, fixing the 
other parameters as in the standard case. 

Figure 3 shows the total fluxes of the forward and reverse shock emission at the observed 
frequency z/ = 5 GHz and the observed times, t = 1 hr, 1 day, 10 days and 100 days, as a 
function of redshift z, together with the 5(7 sensitivities of the VLA and SKA. By comparing 
Figures 3 and 1 we see that the detection is easier than in the standard case, and even the 
VLA could easily detect such energetic GRB afterglows beyond z ~ 30. 

Figure 4 shows the light curves of the energetic GRB afterglows at z = 6 and v ~ 200 
MHz, and at z = 13 and v ~ 100 MHz. The Scr sensitivities of the VLA, LOFAR and SKA 
are also shown. This indicates that LOFAR and SKA can detect such energetic afterglows 
at low frequencies ~ 100 MHz and high redshift z ~ 10, while VLA may detect them if 
observed for a long enough time. The predicted peak flux is about ~ 10-10^ //Jy at the 
redshifted 21 cm frequency, v — 1420/(1 + 2;) MHz, for 2; > 6. 

2.3. High density external medium case 

Even at high redshifts, the ambient density into which the shell expands could be low, 
n < 1 cm~^ as assumed in the standard case (§ 2.1), since the strong radiation from the 

progenitor star may evacuate a cavity in the surrounding gas (Whalen, Abel, & Norman 
2004; Kitayama ct al. 2004). If this docs not happen, however, the external density 
could be higher at high redshifts, based on hierarchical structure formation scenarios. The 
predicted evolution may be, e.g., n oc (1 + zY for a fixed host galaxy mass (Ciardi & Loeb 
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2000). Since little is known observationally about the ambient density at high redshifts, we 
assume a nominal high density case with n — 10^ cm~^. Other parameters are the same as 
the standard case. 

Figure 5 shows the total fluxes of the forward and reverse shock emission at the observed 
frequency u — 5 GHz and the observed times i = 1 hr, 1 day, 10 days and 100 days, as a 
function of rcdshift z, together with the 5cr sensitivities of the VLA and SKA. By comparing 
Figures 5 and 1 we can see that the detection at this frequency is easier than in the standard 
case, and even the VLA can easily detect such high density GRB afterglows beyond z ~ 30. 

Figure 6 shows the light curves of the high density GRB afterglows ai z = 6 and 
u ~ 200 MHz and at ^ = 13 and u ~ 100 MHz. The 5a sensitivities of the VLA, LOFAR 
and SKA are also shown. This shows that the detection of the high density GRB afterglows 
at low frequencies ~ 100 MHz and high redshift z ~ 10 is difficult. Even the SKA would 
require very long observing times. The flux is dimmer than in the standard case in Figure 2, 
because the self-absorption strongly suppresses the low frequency fluxes in this high density 
case. From Figure 6, we see that the peak flux is about ~ 10~^-1 //Jy at the redshifted 21 
cm frequency, u — 1420/(1 + z) MHz, for z > 6. 

2.4. Long duration GRB case 

The GRB duration at high redshifts could in principle be much longer than at low 
redshifts, since more massive progenitor stars could lead to longer accretion times and longer 
jet feeding times (Fryer, Woosley, & Heger 2001). A long duration may also be due to 
continued fallback of material when the explosion is too mild (MacFadyen, Woosley, & 
Heger 2001). BATSE did detect a few bursts longer than 1000 s, with its standard peak- 
fiux trigger, instead of a fiuence trigger. Swift has the ability to operate with a fluence 
trigger, which potentially could lead to the detection of longer GRBs. As an example of long 
duration GRBs, we investigate a case whose duration in the source frame is T/ (1 + ^) = 1000 
s. Other parameters are the same as in the standard case. 

Figure 7 shows the total fluxes of the forward and reverse shock emission at an observed 
frequency z/ = 5 GHz and at observed times, t = 1 hr, 1 day, 10 days and 100 days, as a 
function of redshift z, together with the 5a sensitivities of the VLA and SKA. Comparing 
Figures 7 and 1 we see that the detection of long bursts is easier than for the standard case 
of duration T/ {1 + z) = 10 s, and even the VLA can easily detect such long burst afterglows 
beyond z ~ 30. This is because the reverse shock emission at early times t < 10 day is 
enhanced due to the long duration. 
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Figure 8 shows the long duration GRB afterglow hght curves at 2; = 6 and u ~ 200 
MHz and at ^ = 13 and u - 100 MHz. The 5a sensitivities of the VLA, LOFAR and SKA 
are also shown. Although the reverse shock emission is brighter than in the standard case in 
Figure 2, it is still below the forward shock emission at low frequencies ~ 100 MHz. Since the 
forward shock emission is similar to the standard case, the dctcctability at these frequencies 
is similar to the standard case. LOFAR and SKA should be able to detect the long duration 
GRB afterglows at low frequencies ~ 100 MHz and high redshift z ~ 10, while it appears 
difficult for the VLA to detect them. As seen in Figure 8, the peak flux is about ~ 1-10 //Jy 
at the redshifted 21 cm frequency, u = 1420/(1 + z) MHz, for ^ > 6. 

2.5. Magnetized GRB ejecta case 

Several observations suggest that the ejected shell may be strongly magnetized. A 
large linear polarization in the gamma-ray emission of GRB 021206 has been interpreted as 
implying a uniform ordered magnetic field over the visible region (Coburn & Boggs 2003), 
which could be a magnetic field advected from the central source (see however Rutledge & 
Fox 2004; Wigger et al. 2004). The early afterglows of GRB 990123 and GRB 021211 
also may imply that the reverse shock region has a stronger magnetic field than the forward 
shock region (Zhang, Kobayashi, & Meszaros 2003; Fan et al. 2002b). Although further 
studies are necessary to verify these suggestions (e.g., Zhang & Kobayashi 2004; Matsumiya 
& loka 2003; Fan, Wei, & Wang 2004), it is useful to examine as a possibility the strongly 
magnetized case. Here we adopt a ratio of reverse to forward magnetic field strengths TZb = 5, 
with the other parameters the same as in the standard case. 

Figure 9 shows the total fluxes of the forward and reverse shock emission at the observed 
frequency u = 5 GHz for observed times t = 1 hr, 1 day, 10 days and 100 days, as a function 
of redshift z, together with the 5a sensitivities of the VLA and SKA. Comparing Figures 9 
and 1 we see that the detection of the magnetized GRB afterglow is somewhat easier than for 
the standard case, and even the current VLA could detect such magnetized GRB afterglows 
up to 2; ~ 30. A small increase in the flux is caused by the reverse shock emission at early 
times t <10 day, due to the higher reverse magnetic fleld. 

Figure 10 shows the hght curves of the magnetized GRB afterglows at 2; = 6 and u ~ 200 
MHz and at z = 13 and 1/ ~ 100 MHz. The 5a sensitivities of the VLA, LOFAR and SKA 
are also shown. The reverse shock emission at low frequencies ~ 100 MHz is still below 
the forward shock emission, being dimmer at these frequencies than in the standard case in 
Figure 2. Since the forward shock emission is similar to the standard case, the dctcctability is 
similar to that in the standard case. We see that LOFAR and SKA can detect the magnetized 
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afterglows at low frequencies ~ 100 MHz and high redshift z ~ 10, while the VLA may have 
difficulty detecting them. Prom Figure 10, the peak flux is seen to be about ~ 1-10 //Jy at 
the redshifted 21 cm frequency, v — 1420/(1 + z) MHz, for > 6. 

2.6. Hypernova radio afterglow 

The association of Type Ib/c supernovae (SNe) with GRBs has confirmed the massive 
stellar origin of long GRBs (e.g., Zhang & Meszaros 2003; Meszaros 2002). Only a small 

fraction of core collapse SNe appear to lead to GRB jets with high Lorentz factor 70 > 100 
(Bergcr ct al. 2003; Soderberg, Frail, & Wieringa 2004). However, there may be a larger 
number of SNe accompanied by mildly relativistic ejecta, which could in principle involve 
more energy than typical GRB jets (MacFadyen, Woosley, & Hcgcr 2001; Granot & Ramirez- 
Ruiz 2004). A large kinetic energy ~ 10^^ erg is also suggested in some core collapse SNe, 
by the optical hght curve and spectral fits (Nomoto et al. 2003; Podsiadlowski et al. 2004), 
and these objects appear in a number of cases not to be associated with detected GRB. Such 
energetic supernovae (or hypernovae, HNe) may be powered by jets that fail to bore their 
way through the stellar envelope (Khokhlov et al. 1999; Meszaros& Rees 2001). At high 
redshifts, stars arc expected to be very massive (Abel, Bryan, & Norman 2002; Bromm, 
Coppi, & Larson 2002; Omukai & Palla 2003), and it may be more frequent for the jet to 
be choked even while the jet could be ~ 10 times more energetic (Heger et al. 2003; Fryer, 
Woosley, & Heger 2001). 

Based on these considerations, an isotropic equivalent energy of £^ ~ 10^^ erg may be 
reasonable for high redshift HNe, and depending on the energy conversion efficiency, even an 
isotropic equivalent energy oi E ^ 10^^ erg may not be implausible for population 111 HNe 
(Fryer, Woosley, & Heger 2001; Heger et al. 2003). Keeping in mind the hypothetical nature 
of this assumption, we examine an extreme HNe afterglow case with isotropic equivalent 
energy E — 10^^ erg, ejecta opening half-angle 9 — l/-\/2 (nearly spherical) and initial 
Lorentz factor 70 = 2 (mildly relativistic). Other parameters are the same as in the standard 
GRB case. We note that different authors define HNe differently. Here, our definition of 
HNe is SNe with a large energy and a mildly relativistic ejecta. 

Figure 11 shows the total fiuxes of the forward and reverse shock emission at the observed 
frequency i/ — 5 GHz and observed times t — 1 yr, 3 yrs and 10 yrs, as a function of redshift 
z, together with the 5a sensitivities of the VLA and SKA. This shows that the current VLA 
may be able to detect such energetic HN afterglows up to ^ ~ 20 with ~ 1 day integration 
times, while SKA can easily detect them beyond z ~ 30. We have also calculated a more 
moderate energy HN case with E ~ 10^^ erg, as well as a high external density HN case with 
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n ~ 10^ cm~^. The results for these cases (not shown) indicate that the VLA can detect 
them out to 2; ~ 8 and beyond z ~ 30, respectively, while the SKA can easily detect them 
beyond z ~ 30. 

Figure 12 shows the hght curves of the energetic HN afterglows at 2; = 6 and v ~ 200 
MHz and at z = 13 and v ~ 100 MHz, as well as the 5(7 sensitivities of the VLA, LOFAR 
and SKA. This shows that the VLA, LOFAR and SKA can detect such HN afterglows at low 
frequencies ~ 100 MHz and high redshifts z ~ 10. The peak flux indicated is about ~ 10^- 
10'^ /iJy at the redshiftcd 21 cm frequency, v = 1420/(1 + z) MHz, for z > 6. The reverse 
shock emission is comparable to the forward shock emission, and the rcdshift dependence is 
weak. We also calculated the moderate energy HN case E ~ 10^^ erg and the high density 
HN case n ~ 10^ cm~^ to find that the peak flux at ~ 100 MHz in both cases is about an 
order of magnitude below the flux in Figure 12. 

One caveat about the HN emission is that one may have no prior information about 
the sky position because of the preponderant lack of detected gamma-ray emission. However 
the number of HN afterglows at 2; > 6 on the sky is estimated as ~ 10^ events yr~^ x 
^peak/z/uN/GRB ~ 10^, whcrc tpeak ~ 10^ days is the peak time of the HNe, and we assume 
a fraction ~ 50 % of all GRBs on the sky originating at z > 6 (Bromm & Loeb 2002) 
and a ratio of the HN rate to the GRB one /hn/grb ~ 1- Thus the imaging of ~ 40 deg^ 
of sky would lead to one HN afterglow at z > 6, while the field of view of the SKA may be 
larger than this. Thus, HN afterglows could be detected as a by-product of other related or 
unrelated observing programs. 

It is also interesting to note that the HN afterglows at high redshift may be brighter 
than their low redshift counterparts, since the energy may be larger at high redshifts and 
the redshift dependence of the fiux is weak when the other parameters are fixed. Thus 
the orphan afterglow searches (e.g., Totani & Panaitescu 2002) may preferentially find the 
high redshift afterglows, although the discrimination from variable AGNs may be difficult 
(Levinson et al. 2002). 

3. 21 cm absorption at 2; > 6 

At redshifts z > 6 one may expect 21 cm absorption of the radio afterglow emission 
by the diffuse neutral component of the intergalactic medium (IGM) (Carilli, Gnedin, & 
Owen 2002), as well as by collapsed gas clouds, such as in minihalos and galactic disks 
(Furlanetto & Loeb 2002). The 21 cm absorption by the IGM would provide constraints 
on the reionization history of the universe (Carilli, Gnedin, & Owen 2002), while the 21 cm 
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forest due to gas clouds would provide information about the structure and galaxy formation 
at high redshifts. The 21 cm absorption line due to the host galaxy may provide GRB and 
host galaxy redshifts (Furlanetto & Loeb 2002), independently of other methods such as the 
Lya cutoff (Lamb & Reichart 2000) and X-ray lines (Meszaros& Rees 2003; Gou, Meszaros, 
& Kallman 2004). 

We consider first the optical depth of the diffuse IGM to the 21 cm absorption, which 
is given by 

3cV;i..4i(,//iii(--) n nin^CMB A + ^V^^ ro\ 
^ =^T^-7-277r-E777V = 0-010^^ -TT^ ^H:, (2) 



where Tcmb = 2.73(1+2;) K is the CMB temperature at redshift z, Ts is the spin temperature 
of the IGM, is the neutral hydrogen fraction, /ip is Plank's constant, k is Boltzmann's 
constant, A-^q — 2.85 x 10~^^ s~^ is the spontaneous emission coefficient of the 21 cm line, 
vq — 1420.4 MHz is the source frame 21 cm frequency, H{z) ^ H^VtH^ z)^/'^ is the Hubble 
parameter at high is the neutral hydrogen density at z^ and we assume the helium 

mass fraction to be 24% (Field 1959a; Madau, Meiksin, & Rees 1997). 

The main uncertainty in equation (2) comes from the hydrogen spin temperature T5. 
The spin temperature is determined by the gas kinetic temperature and the coupling of 
Ts with Tk and Tcmb (Madau, Meiksin, & Rees 1997; Tozzi et al. 2000). In the presence 
of only the CMB, we have Tg ~ Tcmb- Once the early generation of stars starts emitting 
Lya photons, the Lya pumping may tightly couple the spin temperature with the gas kinetic 
temperature well before the epoch of full reionization (e.g., Ciardi & Madau 2003). The 
kinetic temperature drops below the CMB temperature without any heating because of the 
adiabatic cooling due to the expansion of the universe, while the kinetic temperature may 
be quickly increased above the CMB temperature by the Lya and X-ray heating due to 
the first stars and the shock heating in the high density regions (Carilli, Gncdin, & Owen 
2002; Chen & Miralda-Escude 2004). The uncertainty in the spin temperature is at least 
0.5TcMB <Ts< 4TcMB (Chen & Miralda-Escude 2004). If Ts ~ 4Tcmb, the optical depth 
is about T ~ 0.0021a;Hi aX. z Q and r ~ 0.0030a;Hi at 2; ~ 13 from equation (2). 

In order to detect the 21 cm absorption hue, the source flux should be larger than 

rmin _ 2A;Tsys SNR 



>leffV2iAi/ r 

where we assume a band width Ai/ ~ 1 MHz, an integration time tint ~ 10 days, the SNR 5, 
the capability of the SKA Aeff/Tgys ~ 5 x 10^ cm^ at i/ ~ 200 MHz, and the smaU r hmit. 
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Therefore ~ 1 mJy sources are needed for the detection of the 21 cm absorption due to the 
diffuse IGM if Ts ~ 4Tcmb- If there was an epoch with low spin temperature Ts ~ O.STcmb 
(Chen & Miralda-Escude 2004), the required flux could be as low as ~ 0.1 mJy, since 
T oc Tg^ from equation (2). Note that the redshift is determined with the precision, 

A..Mi±£)!. 0.070 fi±i)7^y (4) 
I/O V 10 y \lMRzJ ' ^ ^ 

for the band width Au ~ 1 MHz. Note also that the LOFAR will have A^s/Tsys ~ 4 x 10^ 
cm^ at z/ ~ 100 MHz, so that the required flux is ~ 10 times larger than that for the 
SKA. 

Comparing the minimum flux _p™™ of equation (3) with the peak fluxes of the afterglows 
in Figures 2, 4, 6, 8, 10 and 12, we see that in almost all cases there is little or no chance 
to detect the 21 cm line due to the diffuse IGM. The only one exception is the hypothetical 
energetic HN case with isotropic equivalent energy £^ ~ 10^^ erg in a low ambient density 
n < 1 cm~^, if the HN occurs at an epoch with low spin temperature Ts < O.STcmb- 
Although such a HN is not implausible, there is so far no evidence for it. We conclude that 
it is likely to be difficult to detect 21 cm absorption due to the diffuse IGM in the afterglows 
discussed here with the currently existing or envisaged telescopes. Such detections may 
require a telescope with ~ 10 — 100 times better sensitivity than the SKA. 

The 21 cm absorption due to collapsed gas clouds may have better detection prospects, 

since the optical depth can be high, r > 1, because of the high density contrast. The number 
of clouds with optical depth larger than unity t > 1 may be 0.01-0.1 per redshift, while some 
of host galaxies may have r > 1 (Furlanetto & Loeb 2002). In order to resolve the clouds, 
the bandwidth should be narrower than the line width, which is about 

if the line width is determined by the spin temperature. We note that a spin temperature 
of order T5 ~ 10^ K is measured in several damped Lya systems at low redshift, while it 
may be higher at high redshift (Kanekar & Chengalur 2003). If we adopt Au ~ 10 kHz and 
r > 1, the required minimum source flux is about ~ 0.02 mJy from equation (3). Then the 
21 cm absorption line due to gas clouds may be detected in the energetic HN case, as we can 
see from Figure 12, while from Figures 2, 4, 6, 8 and 10, it could be marginally detected in 
the other afterglow cases if the ambient density is low n < 1 cm~^. 
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4. 21 cm absorption at 2; < 6 

At rcdshifts 2; < 6, 21 cm absorption is expected predominantly from neutral gas clouds 
observed as damped Lja systems, rather than from the diffuse IGM, since the latter is 
almost completely ionized. The optical depth of damped Lja systems with the comoving Hi 
column density A^Hi and the spin temperature T5 is about 



3c^hpA,oNn, 

T = 

32'Kki'oTsAi'c 



where we assumed a comoving line width A^c = (1 + z)/^i> — 39(T5/10^K)^/^ kHz as in 
equation (5) (Field 1959a). Thus the optical depth is larger than unity, r > 1, for an Hi 
column density Nn, > 3xlO'^'^{Ts/10^K)^/'^ cm^^. The spin temperature of order Ts ~ 10^ K 
is motivated again on the damped Lya system data and discussion of Kanekar & Chengalur 
(2003). 

As in the previous section, if we adopt Ai/ ~ 10 kHz and r > 1, the required source 
flux is about ~ 0.02 mJy from equation (3). Thus, the 21 cm line due to the damped Lya 
systems may be detected in the energetic HN case, as we see from Figure 12, while it could 
be marginally detected in the other afterglow cases if the ambient density is low n < 1 cm^^, 
from Figures 2, 4, 6, 8 and 10. Note that the detection of the 21 cm line at 2; < 6 is less 
difficult than at ^ > 6, since the observed 21 cm frequency is higher at lower redshifts, and 
the self-absorbed spectrum of the afterglows has a steep rising form, F,, oc u'^ or oc i/^/^, 
in equation (A20). 

The number of absorption lines per unit redshift and unit comoving Hi density is esti- 
mated from observations of the Lya systems as 

where Chj ~ lO^-^ and A ~ 1.33 (e.g., Penton, Stocke, & ShuU 2004), and the redshift evolu- 
tion does not seem to be strong (e.g., Storrie-Lombardi &: Wolfe 2000). Therefore the number 
of the 21 cm line with optical depth larger than unity r > 1 is about ~ 0.006(T5/10^K)^(^~'^)/^ 
per redshift. Since an integration time of order tint ~ 10 days is necessary even for r ~ 1 
clouds, we conclude that it may be difficult to detect the 21 cm line due to the damped Lya 
clouds in the z <Q range. 
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5. Conclusions 

We have investigated the radio afterglows of GRBs and HNe at high redshifts, analyzing 
their detectability as well as their possible usefulness for 21 cm absorption measurements. 
We explored several possible sets of model parameters which are plausible at high redshifts 
(see Table 1). Our source models include the effects associated with reverse as well as 
forward shock emission, jet breaks and sideways expansion of the jet, and the transition 
from a relativistic to non-relativistic expansion regime. 

Our results indicate that standard GRB radio afterglows may be detected with the 
current VLA up to 2; ~ 30, and energetic HN afterglows may be detected up to 2; ~ 20, at 
high radio frequencies ~ 5 GHz. The proposed SKA should be easily able to detect the GRB 
and HN afterglows beyond 2; ~ 30 at frequencies ~ 5 GHz. Other possible effects which may 
be expected at high redshift, such as a larger isotropic equivalent energy, higher external 
density, longer duration GRB and the possibility of magnetized ejecta, can all contribute to 
enhance the detectability of afterglows at ~ 5 GHz. 

We find that it is difficult to detect the 21 cm absorption line due to a diffuse neutral 
IGM, even with the SKA. A possible exception may be in the case of a very high isotropic 
equivalent HN energy E ~ 10"'"^'^ erg, if it occurs in a low ambient density n < 1 cm~^ at an 
epoch with low spin temperature T5 < O.STcmb- This case is very hypothetical, but could 
be plausible under some conditions. 

The 21 cm absorption from collapsed gas clouds along the line of sight, e.g., from 
minihalos and galactic disks at z > 6 and damped Lya systems at 2; < 6, is in principle 
detectable in the energetic HN case, and marginally detectable in the GRB afterglow cases 
considered, if the 21 cm optical depth of the clouds is larger than unity and the ambient 
gas density of the GRBs is low, n < 1 cm~^. However, such high optical depth gas clouds 
are rare, and their detection is not practical given the long integration times needed. This 
indicates that if 21 cm absorption is detected in a narrow range of frequencies, it is likely to 
originate in the host galaxy of the GRB or HN with a high 21 cm optical depth, e.g., in an 
edge-on disk, star-forming galaxy or protogalaxy, and may thus provide a direct measurement 
of the redshift of the host. 

We thank M.J. Rees and S. Kobayashi for useful discussions. This work was supported 
in part by the Eberly Research Funds of Penn State and by the Center for Gravitational 
Wave Physics under grants PHY-01-14375 (KI), NASA NAG5-13286, NSF AST 0098416 and 
the Monell Foundation. 
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A. Afterglow model 

We summarize the afterglow model in this section (c.f., Zhang & Meszaros 2003; Meszaros 

2002) . § A.l is for the forward shock, while § A. 2 is for the reverse shock. A very early evo- 
lution of the reverse shock emission is derived in § A. 3. We will neglect the inverse Compton 
coohng since almost all phase is in the slow cooling (Sari & Esin 2001). We will also neglect 
the effects of the dispersion delay, the free-free absorption, and the scintillation (e.g., loka 

2003) . 

The afterglow evolution is divided into two classes depending on the shell width cT(l + 
z)~^ (Sari & Piran 1995). We consider that the shell is thin (thick) if T < {t^ < T) where 

is the observed time for the forward shock to collect ~ of the shell's rest mass, E — 
10^^ £'53 erg and 70 = 1 0^70,2- The evolution becomes independent of the shell thickness in a 
self-similar manner after the shock crosses the shell at 



tx = max(t^,T), (A2) 

where the first (second) value in the bracket is for the thin (thick) case (Kobayashi & Sari 
2000). The evolution is also independent of the shell width in the very beginning before the 
time 



ti = min(r. In), 

where the first (second) value in the bracket is for the thin (thick) case. 



(A3) 



l + Z 

He 



3E(1 + z) 
167rnmpC^T7Q 



1/2 



(A4) 



is the time for the reverse shock to be relativistic, and T = 10^T2 corresponds to the spreading 
time of the shell (Sari & Piran 1995). Thus the evolution depends on the shell width in the 
interval L < t < ts, . 



A.l. Forwcird shock 



First let us consider the forward shock. We may divide the problem into two parts: 
dynamics and radiation. 
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Dynamics: In the beginning t < U, the swept mass is too small to affect the system and 
the Lorentz factor is constant, 

7(t) ~ 70, (A5) 

as functions of the observed time t. In the interval ti < t < tx, the Lorentz factor is still 
constant for thin shells, j{t) ~ 70, while the deceleration begins for thick shells, 

7(*) ~ 7o(^/^A^)-'/^ (A6) 

since the reverse shock becomes relativistic (Sari & Piran 1995; Sari 1997). After entering 
the self-similar phase, < t, the hydrodynamics is determined by the adiabatic condition 
E ~ iTT'j'^R^nmpC^ /3 and the relation R ~ 47^ct/(l + z) in equation (A13) as 

7(t) ~ [3E(1 + ^)V2567^r^mpC^^Y^^ (A7) 
irrespective of the shell width. The Lorentz factor drops below 6^^ around the time 

to - ^Rg ~ 5 X 10'EliY_/,'n-'/'{l + z) s, (A8) 

where 

Re ~ {3E0yi7mmpc''y/^ (A9) 

Then the jet starts to expand sideways (Rhoads 1999; Sari, Piran, & Halpern 1999), so 
that the dynamics changes as 

7(i) ~ 0-\t/tg)-'/'. (AlO) 

Finally the shock velocity becomes non-relativistic (Waxman, Kulkarni, & Frail 1998; Frail, 
Waxman, & Kulkarni 2000) around the time 



tNR = j^Re ~ 2 X lO^^;^^ + z) s, (All) 



where jj^fR = (1 — P^r)'^^"^ '^^ adopt the non-relativistic velocity Pnr = 2~-^/^. The 
flow becomes nearly spherical, so that the dynamics is described by the self-similar Sedov 
solution as 

m ~ PNRit/tNR)-^^'. (A12) 

For simplicity we use the shock radius given by 

R{t)r.A[^(t)]\t/{l + z), (A13) 
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for the relativistic case (Sari 1998; Waxman 1997; Panaitescu & Meszaros 1998), while 
R{t) (X P{t)t for the non-relativistic case. 

Radiation: Once the dynamics above determine j{t) and R{t), we may estimate the 
synchrotron spectrum as follows. The synchrotron spectrum from relativistic electrons in a 
power-law distribution usually has four power-law segments with three breaks at the self- 
absorption frequency Uaj, the characteristic synchrotron frequency i^mj, and the cooling 
frequency Ucj (Sari, Piran, & Narayan 1998). If we temporarily neglect the self-absorption, 
the observed flux is given by 

(^/^mj)^/^F. ,max,/) ^ < ^mj-i 

'^,1 = S (i^/i^m,/)"^^"^^/^i^z.,max,/, J^mJ < < J^c,f, (AM) 

(^c,//^.n,/)-^-^)/'(^/^c,/)-^/'i^.,ma.,/, ^c,/ < ^, 

in the slow cooling case u^j < Ucj, and 

,max,/) ^ ^ ^c,f^ 

('^/'^c,/)"^/^-^i.,max,/, J^cJ < < l^mj, (Al5) 

,max,/) ^m,f ^ ^) 

in the fast cooling case Ucj < i^mji where 

F„,„„,,= ^^!^^^0i±i) (A16) 

is the observed peak flux at the luminosity distance D — (1 -|- 2;) /q^(1 + z)\dt / dz\dz with 
\dt/dz\-^ = (1 + z)Ho[VLm{l + zf + 0A]^/^ and 

Bf = [327rmpeBjn7(7 - l)]'/'c, (A17) 

is the comoving magnetic fleld. The characteristic synchrotron frequency i/mj and the cooling 
frequency Ucj are given by v{;')mj) and i^ijcj), respectively, where 

= ^f— — , (A18) 

27rmec(l + z)' ^ ' 

is the frequency at which electrons with the Lorentz factor 7e radiate, and 

p-2mp 67rmec(l + ^) 

W = ee,/^-(7-l), 7c,/- ■ (A19) 

Note that we may use Bf and ^rnj in equations (A17) and (A19) even in the non-relativistic 
phase. 



The synchrotron self-absorption hmits the flux below the blackbody emission with the 
electron temperature (e.g., Sari & Piran 1999). This is given by 

F,,BB = 2Tiu^^^,m,{RjD)\l + z)\ (A20) 

where v is the observed frequency, is the observed size of the afterglow, and 7e is the typi- 
cal Lorentz factor of the electrons emitting at v. For simplicity wc use R±{t) — 4'y{t)ct/{l+z) 
for the relativistic case, while R±{t) oc P{t)t for the non-relativistic case. 

A. 2. Reverse shock 

Dynamics: In the initial stage t < ti, the Lorentz factor is constant 7 ~ 70. In the 
interval < t < tx, we have 7 ~ 70 for thin shells, while 7 ~ 7o(^/^iv)~''^^^ for thick shells. 
These evolutions are the same as that of the forward shock. After the shell crossing tx < t, 
the evolution is approximately given by 

7~7x(tAx)^^/^ (A21) 
where 7x = min(7o, Jt) is the Lorentz factor at the shell crossing and 

7t - 100E53/'n-^/^r2"'/'(l + zf, (A22) 

is given by 7(T) from equation (A7) (Zhang, Kobayashi, & Meszaros 2003; Kobayashi 2000). 
A reliable calculation after 7 drops below at te < t requires full numerical simulations 
and is beyond the scope of this paper. Here we use equation (A21) till the velocity becomes 
non-relativistic 7 ~ 7Ar_R, and use the time dependence of the Sedov solution in equation 
(A12) after that. Since the evolution in equation (A21) is the same as that of the sideways 
expansion in equation (AlO), this treatment might not be completely wrong. For simplicity 
we use the shock radius in equation (A13) for all cases. 

Radiation: At the shell crossing t ty^, the forward and reverse shocks have the same 
Lorentz factor 7x, energy density e, and total energy eV. But these shocks have the different 
relative Lorentz factor across the shock 7, i.e., the forward shock has 7/ ~ 7x and the reverse 
shock has 7^ ~ 7o/7x- These shocks may also have different magnetization parameters 
^B,f 7^ ^B,r since the ejected shell may carry magnetic fields from the central source (e.g., 
Zhang, Kobayashi, & Meszaros 2003). Then we can show 

7m,r(tx) 1 7c,r(tx) 1 

rs^ — 

7m,/ (tx) 7^ 7c,/ (tx) 



(A23) 
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(A24) 



^TZb, (A25) 



since 7^ oc 7, 7c oc 7 H ^ and B oc ese^/^, where we define 

7x ^ Br f ^B,r 

— , K-B ^ ^ \ 

70 Bf \€B,f 

We can also show 

^m,r{^x) T^B ^c,riix) 1 -fi/,max,r (^x ) 

I^m,/(^x) 7^' i^c,/(^x) F^,max,/(ix) 

since z/j^ oc jBj^, oc 7^^i3~^i~^, Fj^max c< 'jBNg and A^e e^7~^; where A^e is the total 
number of emitting electrons. 

As in the forward shock case, the synchrotron spectrum is given by equations (A14), 
(A15) and (A20) replacing the subscript / with r, although the emission above u^r disappears 
after the shell crossing < t (see below). Thus the light curves can be calculated by 
specifying the temporal evolutions of Vm, a-nd F,^,max- The evolutions of 7^ and 7c are also 
useful for estimating the self-absorption. In the initial stage t < ti, the evolutions are given 

by 

^m,r OC t , l^c,r ^ ^ ; -P^u,ma,x,r 0^ ^ ) (A26) 
-fm,r OC 7c,r- OC (A27) 

as shown in § A. 3. In the interval ti < t < tx, thin shells have the dependence 

^m,r OC t , l^c,r OC t , -^i/,max,r CC t ^ , (A28) 
lm,r OC t^, -fc,r OC , (A29) 

while the thick shells have 

^m,r OC t , lyf- f OC t , -Pl/,max,r OC t ^ , (A30) 

7^,, <x 7c,, oc t-^'\ (A31) 

After the shell crossing < t, no electron is shocked and electrons only cool adiabatically, 
so that the evolutions are given by 

^m,r OC i , Vc,r OC t , -^i/,max,r OC t , (A32) 
-1m,r OC 7c,. oc r^'\ (A33) 

and the emission above z/c,, disappears (Zhang, Kobayashi, & Meszaros 2003; Kobayashi 
2000; Gou, Meszaros, Abel, & Zhang 2004). In the Sedov phase we use 

V^^r OC t~^, Vc,r OC Fj.,max,r OC fr^''\ (A34) 

7„,,ocr^/^ 7c,.oc^-^/^ (A35) 
by noting that the energy density drops as oc oc t~^/^. 



-20- 



A.3. Very early reverse shock 

Let us derive the evolution of the reverse shock emission in the initial stage t < tj. The 
evolution is independent of the shell width. The Lorentz factor is constant 7 ~ 70 since 
the reverse shock is non-relativistic 7^ — 1 <^ 1 (Sari & Piran 1995). This together with 
the constant ambient density leads to a constant energy density in the forward and reverse 
shock Cr ~ const. Since the shell spreading has not started, the shell density decreases as 
rir oc oc t'"^. Then the relative velocity across the reverse shock is Pr (7^ — ly^'^ oc 
ey^rir oc t, and hence the total number of emitting electrons increases as A^e oc f3rt oc t^. 
Recalling 7^ oc Cr/rir, 7c oc ^~^B~H~^ and B oc 6^^, we have the time dependence in 
equation (A27). Similarly we can derive equation (A26) since oc jBj^, Uc oc ^y^^B'H'^ 
and F^,max oc -fBNg. 
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Table 1. Model parameters and some results: the isotropic equivalent energy E erg, the 

opening half-angle 6, the ISM density n cm~^, the spectral index p, the initial Lorentz 
factor 7o, the duration in the source frame T/ (1 + z) s, the energy fraction of electrons eg, 
the energy fraction of magnetic fields e^, the ratio of reverse to forward magnetic field TZb, 
the maximum redshift out to which these afterglows could be detected with the VLA zvla, 
and the peak flux at the redshifted 21 cm frequency, u — 1420/(1 + z) MHz, for z>6, 

1+z 
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E [erg] 
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70 
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Fig. 1. — GRB afterglow total flux from the forward and reverse shock emission at an 
observed frequency z/ = 5 GHz and observed times 1 hr, 1 day, 10 days and 100 days, as a 
function of redshift z. The standard model parameters used are: isotropic equivalent energy 
E = 10^'^ erg, opening half-angle 6 = 0.1, ISM density n = 0.1 cm~^, electron spectral index 
p — 2.2, initial Lorentz factor 70 = 200, GRB duration in the source frame T/(l + z) ~ 10 
s, and plasma parameters €e — 0.1, = 0.01, TZb — 1- The 5a sensitivities of the VLA 
(dotted line) and SKA (dashed line) for an integration time tint — 1 day and a band width 
Au — 50 MHz are also shown. 
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Fig. 2. — GRB afterglow forward shock (solid line) and reverse shock (dashed line) fluxes, 
shown as a function of the observed time t at frequencies near the redshifted 21 cm radiation 
for two representative redshifts. We used the standard GRB model parameters: isotropic 
equivalent energy E = 10^^ erg, opening half-angle 6 = 0.1, ISM density n = 0.1 cm~^, 
electron spectral index p = 2.2, initial Lorentz factor 70 = 200, GRB duration in the source 
frame T/{1 + z) — 10 s, and plasma parameters eg = 0.1, = 0.01, TZb — 1- The thin hnes 
are for a redshift z — 6 and observed frequency 1^ — 200 MHz, while thick hnes are for 2; = 13 
and v — 100 MHz. The 5(7 sensitivities of the VLA, LOFAR and SKA for integration times 
of one-third of t and band width Au — 50 MHz are also shown by the long dashed lines. 
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Fig. 3. — GRB afterglow fluxes, same as in Figure 1 except for the isotropic equivalent 
energy, which is taken here 0,8 E — 10^^ erg. 
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Fig. 4. — GRB afterglow fluxes at two redsliiftcd 21 cm frequencies, as in Figure 2 except 
for the isotropic equivalent energy, which is taken here as = 10^^ erg. 
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Fig. 5. — GRB afterglow fluxes as in Figure 1 except for the ISM density, which is taken 
here as n = 10^ cm~^. 
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Fig. 6. — GRB afterglow fluxes at two redshifted 21 cm frequencies, as in Figure 2 except 
for the ISM density, which is taken here as n = 10^ cm~^. 
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Fig. 7. — GRB afterglow fluxes as in Figure 1 except for the GRB duration in the source 
frame, which is taken here as T/ z) — 10^ s. 
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Fig. 8. — GRB afterglow fluxes at two redshifted 21 cm frequencies, same as in Figure 2 
except for the GRB duration in the source frame, which is taken here as T/ {1 + z) ~ 10^ s. 
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E=io53erg, = 0.1, n = 0.1cm-3, p = 2.2 
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Fig. 9. — GRB afterglow fluxes, as in Figure 1 except for the magnetic fleld in the reverse 
shock, which is taken here as TZb — if-B,rl ^bjY^"^ — 5- 
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Fig. 10. — GRB afterglow fluxes at two redshifted 21 cm frequencies, as in Figure 2 except 
for the magnetic field in the reverse shock, which is taken here as TZb — {^B,r/^B,fY^'^ — 5. 
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Fig. 11. — Hypernova (HN) remnant total flux from the forward and reverse shock at the 
observed frequency i/ — 5 GHz and observed times 1 yr, 3 yrs and 10 yrs, shown as a function 
of redshift z. The model parameters are as for the standard GRB case in Figure 1, except 
for the isotropic equivalent energy E = 10^'' erg, the opening half-angle 9 = 2 (nearly 
spherical) and the initial Lorentz factor 70 = 2 (mildly relativistic). The 5a sensitivities of 
the VLA (dotted line) and SKA (dashed line) for an integration time tint = 1 day and a 
band width Au = 50 MHz are also shown. 
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Fig. 12. — Hypernova (HN) remnant fluxes at two redshifted 21 cm frequencies, same as 
Figure 2 except for the isotropic equivalent energy E — 10^^ erg, the opening half-angle 
9 — 2^^/^ (nearly spherical) and the initial Lorentz factor 70 = 2 (mildly relativistic). The 
5o" sensitivities of the VLA, LOFAR and SKA for the integration time of one-third of t and 
the band width Ai/ = 50 MHz are also shown by the long dashed lines. 



